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ABSTRACT 

We have developed a simulation model of particle acceleration in coronal shock waves. The 
model is based on a Monte Carlo method, where particles are traced in prescribed large-scale 
electromagnetic fields utilizing the guiding center approximation. The particles are scattered in 
the turbulence according to quasilinear theory, with the scattering amplitude directly proportional 
to the intensity of Alfven waves at gyro-resonant wavenumbers. The Alfven waves are traced 
simultaneously with the particles, so that the wave field is propagated outwards from the Sun 
using WKB propagation supplemented with a phenomenological wavenumber diffusion term and 
a growth rate computed from the net flux of the accelerated particles. We consider initial wave 
amplitudes small enough to allow rapid escape of particles from the shock to the ambient medium. 
Thus, in our model the Alfven waves responsible for the diffusive acceleration of particles are 
generated by the accelerated particles themselves. In this work, we study the effects of non- 
constant shock velocity and non-monotonic Alfven velocity on particle acceleration scenarios. 
We report in particular how the deceleration of a shock affects particle intensity and turbulence 
power evolution in the vicinity of the shock. 

Subject headings: coronal shock particle acceleration, decelerating shock 
PACS: 96.50.Pw 96.50.Tf 96.60.ph 



1. Introduction 

Acceleration by coronal and interplanetary 
shocks driven by coronal mass ejections (CMEs) 
is currently considered the primary source of large 
solar energetic particle (SEP) intensities. Injected 
particles repeatedly cross the shock front, scat- 
tering from the plasma waves, and gain energy 
on each crossing. A requirement for rapid ac- 
celeration and escape is strong turbulent trap- 
ping in front of the shock and weaker trapping 
further out. Such a turbulent structure can, 
in large events, be maintained by streaming of 
the accelerated particles themselves, as plasma 
waves in the upstream experience net amplifi- 
cation by scattering events. Diffusive accelera- 
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tion of solar energetic particles in coronal shocks 
driven by fast CM E s has been discussed by eg . 



Axford et al. (1977); Blandf ord fc Ostrikerl ( 197 



and lKrvmskiil ([19771 ). wher eas upstrea m wave am- 
plification was detailed by iBeTi (119781) . 

The focus of recent research has been timc- 



depen dent numerical modeling. IVainio fc Laitinen 
(|2007l ) reported on an approach using a simplified 
solar wind model utilizing constant wave group 
speed V — u sw + va and constant shock veloc- 
ity. The initial study, while bringing new insights 
into the energetic particle research, could not di- 
rectly be extended to variable Alfvenic speeds as 
the model used a Lagrangian grid with tempo- 
rally constant spacing to simplify wave transport 
modelling. 

In this study, we use a shock-attached wave 
grid, which allows realistic coronal velocity pro- 
files and greater numerical accuracy near the 
shock. In addition, eliminating interpolation cal- 
culations at the shock front stop numerical oscil- 
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lations from forming. Our model will allow us 
to implement multiple (in- and outbound) wave 
modes in future studies. Spatial convection of 
wave power within the grid was done wi th a Van 
Leer flux-limited Lax-Wendroff scheme (van Leer 



1974 : lLax fe Wendr off 1960). 



2. model 

2.1. Wave-particle interactions 

In our coupled SEP accele ration and Alfven 
wave self-generation model ([Vainio fe Laitinen 
20071 ). the particles are simulated with a Monte- 
Carlo method, while the waves are described with 
a frequency power spectrum, propagated in the 
plasma background with the WKB wave trans- 
port equation supplemented by spectral transport 
and amplification terms, 



dP dP ~ 

h v = aP ■ 

dt dr 



We use a simplified form of the amplification co- 
efficient 



2 , pS P {r,t) 



uva 



as derived in IVainid (12003b . where 

S p (r,t) = 2ttp 2 v fiF(r,p 7 fi,t)dfi 



represents the particle streaming in the Alfven 
wave frame. For the wave diffusion coefficient, 
we use a linear diffus ion model as given in 
Vainio fe Laitinen! (|2007l ). 



Dff = (V/r e )f/ 3 f; 2/3 . 

In the above equations, P is the normalized wave 
power P = (V 2 /Bva)P, f cp is the proton cy- 
clotron frequency and p = rn p Vf cp /f is the reso- 
nant wave- frame momentum simplified by neglect- 
ing its dependence on particle pitch angle. The 
wave group speed is V, r© = 1 AU, F is the 
proton distribution function, and the breakpoint 
frequency fb is set to 1 mHz to match detected 
spectra at 1 AU. Streaming is tracked via particle 
movement between simulation cells, weighted to 
take into account the movement of the grid. Wave- 
particle interactions are implemented so that par- 
ticle energies can be relativistic. Particles are not 



traced downstream of the shock, but rather in- 
stantly returned upstream with their statistical 
weight adjusted by the theoretical return proba- 
bility. 

2.2. Solar wind model 

Our solar wind model consists of a back- 
ground solar wind with speed u(r) 7 electron den- 
sity n(r) and a super-radially expanding coro- 
nal/interplanetary magnetic flux tube. The mag- 
netic field at the axis of the flux tube is 



*(•-> = *(t)'" 



1 + 1.9I& 

r 



with Bo scaling th e values to 2.90 nT at 1 AU 
( Vainio et al.ll2003h . The cross-section of the flux 
tube is inversely proportional to the magnetic 
field, and the off-axis component of the mag- 
netic field, while satisfying the V • B = con- 
dition everywhere, vanishes on the axis. The 
parallel shock is taken to propagate at speed 
U s (r) along the mean magnetic field. For the 
electron density, we use the empi rical model in 
Cranmer fe van Ballegooiienl ([20051 ) . 



n(r n ) = n (1.0 r- 2 + 25.0 r,; 4 + 300.0 r,; 8 
+1500.0 r- 16 + 5796.0 r" 33 9 ), 

with r n = r / Rq and n set to facilitate, with the 
rest of the solar wind model, an Alfven speed of 20 
km/s at 1 AU. The solar wind velocity is inferred 
from mass conservation. This solar wind model 
results in a local maximum in the Alfven velocity 
below 5Rq. 

2.3. Simulation parameters 

In this study, we propagate shocks from 1.5 Rq 
with an initial velocity of 1500 km/s. Turbulence 
is tracked on a logarithmic grid reaching out to 
300 Rq in front of the shock. The mean magnetic 
field is normal to the plane of the shock, and the 
scattering center compression ratio is set to 4. The 
ambient 100 keV proton mean free path is set to 
Ao = IRq at ro = 1.5Rq. The total amount of 
injected particles (within a flux tube of 1 sr solid 
angle at the solar surface) is given values of 2 • 10 35 
(run A), 6 • 10 35 (run B) and 2 • 10 36 (run C). The 
injection of a sligh tly suprathermal particle popu- 
lation is detailed in IVainio fc Laitinenl (|2007l ) . The 
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Fig. 1. — Solar wind and shock parameters; the 
line for the constant velocity shock speed is trun- 
cated so it does not overlap with the legend. 

shocks are propagated with both constant velocity 
and with a constant deceleration of 50 m/s 2 . 

3. Results and discussion 

Our simulations show that both constant veloc- 
ity and decelerating shocks are capable of signif- 
icant acceleration of energetic particles, although 
constant velocity shocks do attain greater maxi- 
mum energies. For run A, the gradually decreas- 
ing shock velocity proves to be significantly less 
effective in accelerating particles. For run B, the 
decelerating scenario retains greater turbulence in 
front of the shock, which results in stronger trap- 
ping and continued acceleration. This simulation 
displays a continuous proton power-law spectrum 
up to 1 MeV, whereas the constant velocity case 
retains less high-energy particles near the shock. 
The increased acceleration efficiency and resulting 
turbulence generation of case C causes the trap- 
ping difference to become insignificant, with both 
constant velocity and decelerating shocks display- 
ing a continuous proton power-law spectrum up to 
20 MeV. 

3.1. Comparison at set propagation dis- 
tance 

The particle and wave turbulence distributions 
(figures [3] and 02 respectively) resulting from run 
B are presented when both shocks have travelled 
the set distance of 22.7 R®. 
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Fig. 2. — Proton intensities with injection of 
6 • 10 35 particles (run B), when r s = 22.7 R Q . 
The constant velocity shock values are represented 
with contour lines, whereas the decelerating case 
values are presented with filled contours. Contour 
lines are at one magnitude intervals. 

Particles from the decelerating shock have trav- 
elled further due to greater total time simulated, 
but have not reached as high energies as those ac- 
celerated by a constant velocity shock. The en- 
hancement in turbulence in front of the decelerat- 
ing shock is stronger and more diffuse than that 
of the constant velocity shock, which can be ex- 
plained by prolonged exposure to efficient wave 
generation at lower altitudes, and slower convec- 
tion of the upstream plasma to the shock. 

3.2. Comparison at set propagation time 

The particle and wave turbulence distributions 
(figures |4] and [5l respectively) resulting from run 
B are presented when both shocks have travelled 
the set time of 1 64 minutes. Clearly, the constant 
velocity shock accelerates particles to higher max- 
imum energies and the particles travel further. 

3.3. Spectra in front of the shocks 

Figures |6l [7] and [8] show the energetic parti- 
cle intensities and wave power spectra found in 
front of the shocks for runs A, B and C, respec- 
tively. Each figure displays both the spectrum for 
a constant velocity shock and spectra for deceler- 
ating shocks, when the decelerating shock has ei- 
ther travelled the same time or the same distance 
as the constant velocity shock. For comparison, we 
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Fig. 3. — Wave powers with injection of 6 ■ 10 35 
particles (run B), when r s — 22.7 i?©. For contour 
descriptions, see figEl 



Fig. 5. — Wave powers with injection of 6 • 10 35 
particles (run B), when t — 164 minutes. For con- 
tour descriptions, see figEl 
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Fig. 4. — Proton intensities with injection of 6 • 
10 35 particles (run B), when t = 164 minutes. For 
contour descriptions, see figEl 



Fig. 6. — Proton intensities and wave powers in 
front of the shock, injection of 2 • 10 35 particles 
(run A). 



also present th e spectrum according to the steady- 
state model of Bell as further explained in 
IVainio fc Laitinenl (|2007l ). The upstream plasma 
velocities (927 and 517 km/s) are calculated from 
the solar wind and shock velocities at the highest 
altitudes attained (11.5 R@ and 22.7 Rq, respec- 
tively). It should be noted, that the upstream 
plasma velocity does not remain constant in our 
simulations as they do not represent a steady-state 
scenario. 



4. Conclusions 

Our results suggest that decelerating shocks can 
accelerate particles to high energies, though not 
as efficiently as in the case of a constant velocity 
shock. The acceleration efficiency and the nature 
of particle trapping vary greatly based on the num- 
ber of injected particles, which proves that further 
parameter studies are required to shed light on the 
dynamical, non-linear effects surrounding deceler- 
ating shocks. In addition, the effects of the realis- 
tic evolution of the scattering center compression 
ratio and the particle injection profile based on the 
solar wind density remain an open question. 
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Fig. 7. — Proton intensities and wave powers in 
front of the shock, injection of 6 • 10 35 particles 
(run B). 
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Fig. 8. — Proton intensities and wave powers in 
front of the shock, injection of 2 • 10 36 particles 
(run C). 
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